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1. BACKGROUND 

At present there is no single “best way” to predict propagation 
impairments to an Earth-Space path. However there is an internationally 
accepted way — namely that given in the most recent version of CCIR 
Report 564 of Study Group 5. This paper treats a computer code 
conforming as far as possible to Report 564. It was prepared for an IBM 
PS/2 using a 386 chip and for Macintosh SE or Mac II. It is designed to be 
easy to write and read, easy to modify, fast, have strong graphic 
capability, contain adequate functions, have dialog capability and windows 
capability. 

Computer languages considered included the following 

•Turbo BASIC 
•Turbo PASCAL 
•FORTRAN 
•SMALL TALK 
•C++ 

•MS SPREADSHEET 
•MS Excel-Macro 
•SIMSCRIPT II.5 
•WINGZ 

Microsoft Excel-Macro was chosen as the first phase simulation language 
for the following Characteristics 

• strong graphic capability 

• about 400 math or control functions 

• sophisticated coordinate systems 

• window and dialog capability 
•easy to customize 

• enough resolution and colors for our use 

• not too fast at the beginning, because of on-screen processing 
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2. PROGRAM STRUCTURE 

Shown in Figure 1 are the dialog boxes illustrating how the first part of 
the program would be run. The second part consists of consolidating 
gaseous attenuation and rain attenuation with the free-space value of the 
carrier C; and brightness temperature with the noise temperature of the 
receiving system, to obtain the overall system noise density Xq or Nq to 



dialog boxes 


Figure 1 Program Diagram 
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3. METHOD AND RESULTS 

(A) Rain Rate as a Function of percentage of year: CCIR Report 

563 (1990) is used as the authority for rain rate. The formula used is 

-ur 

P(R>r) = — — r>2 mm/h (1) 

b 
r 

where R is rain rate, r is a given rain-rate threshold, u is a parameter 
depending on climate and geographical features, and 
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Figure 2 illustrates this relation applied to 4 CCIR regions (found in Report 
563), and Figure 3 illustrates the goodness of fit of equation (1) with the 
CCIR data in Report 563. As can be seen, beyond 2 mm/h the fit appears 
very good. 

(B) Attenuation due to Rain: CCIR Report 564 (1990) is used to 

obtain rain attenuation. The specific attenuation Yr is obtained from 

Y r = k (Rq q|) dB/km (2) 

where k and a are tabulated values and are a function of frequency and drop 
size distribution. The path attenuation A is then given by 

A = T r L s r 0.01 dB (3) 


where L s is the adjusted path length through rain, Ro.oi is the rain rate in 
mm/h for 0.01% of the time. Figure 4 illustrates these relations in a plot of 
attenuation vs frequency for four different rain rates. 

(C) Gaseous Attenuation: According to the CCIR Report 719 

(1990) version, specific attenuation due to oxygen and water vapor are 
determined as 
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where y 0 and y w represent specific attenuation by oxygen and water 
vapor, f is frequency in GHz, and p w is water vapor density. Total 
attenuation A g due to gaseous then is determined as 
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Figure 5 illustrates attenuation due to the gaseous atmosphere at four 
different elevation angles, and Figure 6 shows attenuation due to the 
gaseous atmosphere at four different station heights. 

D: Brightness Temperature (upward looking antenna): The 

general formula used for brightness temperature was based on the model 
developed by Waters in 1976. Because it is quite complicated, we used a 
simplified formula which was the combination of two models developed by 
E. K. Smith (1982) and Waters (1974) 


A = 



[7o+Y w l 


dr 



-T(r) , -t„ 

T(r) y(r) e dr+T^e 


where 


T(r) = 
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If T(r) is replaced by a mean path temperature Tm, 


simplified as 


T b =T m <i-e ) 


it can be 


or 


T b =T m (l-10 [ - A(dBV10sin91 ) 


e>io 


O 


(9) 


where A is the attenuation at zenith direction, Tb is brightness 
temperature* 9 represents elevation angle, L is the loss factor, and x is 
the optical depth. 

Figure 7 illustrates the brightness temperature at four different 
elevation angles, and Figure 8 is a comparison of computed results and 
CCIR data at the conditions of water vapor density 7.5 mm/h and elevation 
angle 30 degrees. As we can see, the fit is very good in lower frequency 
range. 
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Figure 2 Rain Rate vs % of Year for Regions B, E, K, M 
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Figure 3 Analytical vs Experimental Results (CCIR region E) 
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Figure 4 Specific Attenuation y r Due to Rain 
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Figure 5 Gaseous Attenuation Transiting the Atmosphere 
(pressure 1013 mb, temperature 15° C- Water Vapor 7.5 g/m3) 
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Attenuation (dB) 
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Brightness Temperature (K) 
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Figure 7 Brightness 







Figure 8 Brightness Temperature Comparison of the Simple Model 
and the CCIR Radiative Transfer Model 
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The program takes about 2-5 minutes to run for each application, 
depending on the parameters given. The program can be totally 
customized or manually controlled. In the former case, the program can 
open, run, and close automatically; the only thing needed to do is to input 
your options and answer questions asked. 
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Software for Propagation 
by W. Vogel, U. of Texas, Austin 


1 . Choice of Hardware and Software Platform 

2. Selection of Problems To Be Coded 

3. Coding & Solicitation of Contributions 

4. Testing & Documentation 

5. Dissemination to Users 

6. Technical Support 

7. Revisions/Expansion 


1. Hardware and Software Platform 


• IBM/PC + Macintosh 

• Spreadsheet: 

• Math: 

• Special Purpose: 

• Language: 

• Survey Users 

• User Modifiable 


1-2-3, EXCEL, QUATTRO... 

MATHCAD 

EE-PAC 

Quick, FORTRAN 
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2. Problem Selection 

• CCIR Greenbook 

• NASA Handbook 

. Books: P+B, Ippolito... 

• Models 

• Orbits 

• Footprints 

• Propagation Database 

• Call for Contributions 


3. Coding & Solicitation of Contributions Standards 

4. Testing & Documentation 
Operation, Logical Errors, B-Testing 
Easy To Use, Built-In Help, Example 
User Knowledgeable 


5. Dissemination 
To Anybody 
Mail, BBS 
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6. Technical Support 
BBS, FAX 

7. Revisions/Expansion 
Useful Lifetime... 

Enhance Functionality 
Include New Developments 
Advisory Panel 
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